Mass flux characteristics in solid 4 He for T > 100 mK: Evidence for Bosonic Luttinger 

Liquid behavior 
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At pressure ~ 25.7 bar the flux, F, carried by solid 4 He for T > 100 mK depends on the net 
chemical potential difference between two reservoirs in series with the solid, Aii, and obeys F ~ 
(A/i) 6 , where b ~ 0.3 is independent of temperature. At fixed A/x the temperature dependence of 

the flux, F, can be adequately represented by F r ln(T/r), r ss 0.6 K, for 0.1 < T < 0.5 K. A 

single function F = Fo(Afi) b \n(T/r) fits all of the available data sets in the range 25.6 - 25.8 bar 
reasonably well. We suggest that the mass flux in solid 4 He for T > 100 mK may have a Luttinger 
liquid-like behavior in this bosonic system. 

PACS numbers: 67.80.-s, 67.80.B-, 67.80. bd, 71.10.Pm 
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Following the measurements of Kim and ChanQ, i 
and the interpretation of the possible existence of a 
supersolid 3] , there has been renewed interest in solid 
4 He. Some have questioned the supersolid interpretation 
and imply that some experiments carried out to date may 
show no clear or only weak direct evidence for supersolid 
behavior 0, Experiments designed to create flow in 
solid 4 He in confined geometries by directly squeezing 
the solid lattice have not been successful [6 -9]. We took 
a different approach and by creation of chemical poten- 
tial differences across bulk solid samples in contact with 
superfluid helium have demonstrated mass transport by 
measurin g th e mass flux, F, through a cell filled with 
solid 4 He T3, U| at temperatures that extend to values 
above those where torsional oscillator or other experi- 
ments have focused attention. Indeed these experiments 
revealed interesting temperature dependence 12|, Il3j in 
the vicinity of 80 mK, where the major changes in tor- 
sional oscillator period or shear modulus fl4| were seen. 

Here we seek to understand the behavior of F for T > 
100 mK in more detail. We apply a temperature dif- 
ference, AT, to create an initial chemical potential dif- 
ference, A/tto, between two superfluid- filled reservoirs in 
series with a cell filled with solid He. We then mea- 
sure in some detail the behavior of the 4 He flux through 
the solid-filled cell for T >100 mK that results from the 
imposed AT as the pressure difference between the two 
reservoirs changes (the fountain effect) and the chemi- 
cal potential difference between the two reservoirs, A/z, 
changes from A/j, to zero. For T > 100 mK modest pe- 
riod shifts have been seen in a number of torsional oscil- 
lator experiments, in some cases even above 400 mK[l5j|. 

Since the apparatus [l2l Fui ] used for this work has been 
described in detail previously, our description here will be 
concise. A temperature gra dient is present across the 
superfluid-filled Vvcor [17H19j rods (Figure 1), VI and 
V2, which ensures that the reservoirs Rl and R2 remain 
filled with superfluid, while the solid-filled cell (1.84 cm 3 ) 
remains at a low temperature. For the present experi- 



ments a chemical potential difference can be imposed by 
the creation of a temperature difference, AT = Tl — T2, 
between the two reservoirs. The resulting change in the 
fountain pressure [2^] between the two reservoirs results in 
a mass flux through the solid-filled cell to restore equilib- 
rium. The experimental protocol is designed to minimize 



what has been described as the "syringe effect" [21|, |2 



by which sequential net injections of atoms to the cell 
increase the density of the solid. By a reduction in the 
base temperatures of Rl and R2 we also eliminate the 
flow restriction that would be present for too high a Vy- 
cor temperature [13j. 

To fill the cell initially, the helium gas (ultra high pu- 
rity; assumed to contain ^300 ppb 3 He) is condensed 




FIG. 1. (color online) Diagram of the apparatus. The pres- 
sure of the solid is measured by capacitance strain gages [IB] 
CI and C2, the pressures in the superfluid-filled reservoirs, 
P1,P2, are measured at room temperature, and the sample 
cell temperature is measured by thermometer TC. [Not to 
scale; VI and V2 are longer than shown here.] 
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through a direct-access heat-sunk capillary (not shown 
in Figure [lj . To grow a solid at constant temperature 
from the superfluid, which is our standard technique, we 
begin with the pressure in the cell just below the bulk 
melting pressure for 4 He at the growth temperature and 
then add atoms simultaneously through lines 1 and 2. 
Once we have created solid at the desired pressure, we 
close the fill lines and change the cell temperature. 

With stable solid 4 He in the cell, we use heaters HI 
(H2) to vary Tl (T2) to create chemical potential differ- 
ences between the reservoirs and then measure the result- 
ing changes 20] in the pressures PI and P2. An exam- 
ple of the behavior seen from an application of this ap- 
proach is shown in Figure EJ where PI, P2, AP = P1 — P2 
and Tl and T2 are shown as a function of time. A 
baseline reservoir temperature is first selected, To, with 
Tl = T2 = T . Then Tl is decreased by ST while T2 is 
increased by the same interval; AT = Tl — T2 = —2ST. 
After chemical potential equilibrium is reached (A, Fig- 
ure 2), the values of Tl and T2 are interchanged (B, Fig- 
ure 2); AT = Tl -T2 = +2ST. Later ST is changed by a 
small amount (C, Figure 2) and the process is continued 
as time elapses. With each switch in the value of Tl — T2 
there is a response of Tl — P2. This approach is expected 
to create a smaller perturbation on the solid and allows 
us to obtain larger Ali values without exceeding the up- 
per Vycor temperature at which a significant flow limita- 
tion is encountered H3- We take F = d(Pl - P2)/dt to 
be proportional to the flux of atoms that passes through 
the solid. We study F as a function of T and All, the 
chemical potential difference between Rl and R2, where 
All = rriilj (dP/ ' P) — / (sdT)], where rri4 is the 4 He mass, 
p is the density and s is the entropy per unit mass. We 
report All in units of J/g instead of J/atom. We will 
report our flux values in mbar/s, where a typical value of 
0.1 mbar/s corresponds to a mass flux through the cell 
of » 4.8 x 10~ 8 g/sec. 

We typically consider the data in two ways: (1) F as 
a function of AT at a sequence of fixed solid 4 He tem- 
peratures and (2) F for fixed AT as a function of T. For 
data of the first sort, we measure the dependence of the 
flux F, on the imposed temperature difference between 
the reservoirs Rl and R2, AT. Following the applica- 
tion of the imposed AT the system responds with a flux 
to create an increasing AP (the fountain effect). Thus 
the net chemical potential difference A/i between the two 
reservoirs decreases to zero as AP increases. Since the 
flux should depend on Ali we document that behavior. 
An example of the relationship between F and All is 
shown in Figure [3] for several solid 4 He temperatures. 
These data have error bars that are related to our ability 
to determine the flux from the measured d(Pl — P2)/dt 
and this becomes more difficult at small values of A/i. 
We find that a reasonable characterization of the data is 
given by F — ^4(A/i) b . The results of fitting the data 
to this functional form for several sets of data at various 
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FIG. 2. (color online) Response of pressures PI and P2 
to the application of two different ST. T = 390 mK. Use of 
heaters, HI and H2, results in changes in Tl and T2. The 
resulting changes in PI and P2 are best seen as PI — P2, 
shown here (uppermost data). The small drift in PI, P2 of 
the sort seen here is typical and variable and appears to have 
no influence on PI — P2 



cell temperatures are shown along with the data in Fig- 
ure [U We find that A has temperature dependence, but 
that the exponent b is constant within our errors, as is 
illustrated in Figure 2] for several data sets. We conclude 
from the behavior seen in Figure [3l which we have seen in 
other samples, that our measurements here are primar- 
ily in the dissipative regime. This dissipation may come 
from phase slippages, a many-body tunneling phenomena 
expected in a superfluid-like system. Historically, some 
have explored the approach to the dissipative regime in a 
superfluid system by study of the flow velocity associated 
with a pressure gradient [23] or a decreasing gravitational 
pressure head[24j[. under quasi- isothermal conditions. An 
interchange of the axes of Figure [3] is reminiscent of such 
studies. 

Next, we study F vs. T for several fixed values of A/i. 
We find that the data of this sort can be reasonably well 
represented by a function of the form F ~ — ln(T/r), 
where r is a fitting parameter. As a specific example, 
in Figure [5] we show data deduced from that shown in 
Figure 3 for which we find that r w 630 ± 20 mK, 
which is consistent with earlier observations in this pres- 
sure range fl3j. which showed no evidence for flux above 
w 650 mK. Typically in the temperature range studied 
here the behavior of the flux is reproducible and not hys- 
teretic; for a given value of A/i one can reproduce the 
measured F(T) value (e.g. shown in Figures 3, [5]) for 
increases or decreases in temperature. But, at times the 
flux becomes unstable and can fall to values which are 
indistinguishable from zero. Once this happens, or the 
temperature is raised above ~ 650 mK and then lowered, 
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FIG. 3. (color online) Values of F = d(Pl - P2)/dt for 
a solid sample at a cell pressure of (CI + C2)/2 = 25.6 - 
25.8 bar shown as a function of Afi determined for the case 
AT = 27 mK; | AT |= 25T. Small shifts in A/x have been 
applied to align the data at A^t = 0. The non-hysteretic flux 
depends on A/i, and can be represented by F — A(A(i) b . 



it almost never is the case that a finite value of the flux 
will reappear when the temperature is lowered unless net 
atoms are added to or removed from the cell, which pre- 
sumably changes the disorder in the solid. 

We find that a single function F = F (A^) f ' ln(T/r) 
fits the data in the range 25.6 - 25.8 bar reasonably well. 
From simultaneous fits to the data for the dependence on 
AyU and T for the data of Figures 3 and 5, we find average 
values for the parameters to be Fq = — 0.50±0.03 mbar/s, 
b = 0.29 ± 0.01 and r = 0.63 i 0.01 K. 

The behavior of F vs. A/i appears to be qualitatively 
different from the behavior seen by Sasaki et al.[2^| for 
flow along grain boundaries on the 4 He melting curve; 
our flux decreases with decreasing A/i, theirs did not. In 
a saturated vapor pressure study of the decay of a super- 
fluid 4 He level, h, due to flow through narrow (~ 2-5 /im) 
slits formed between two flat plates Rorschach found 
that dh/dt = ah 1 / 3 , a dependence on h consistent with 
Gorter-Mellink friction. Presumably A/i ~ h in each 
case. The temperature dependence found by Rorschach 
closely followed the temperature dependence of the bulk 
superfluid density, which is distinctly different from the 
temperature dependence we observe, F ~ — m(T/r). 

Given the fact that we have a solid-filled cell off the 
melting curve, an unusual dependence on T, dissipa- 
tive flux and sensitivity to disorder, we suggest a pos- 
sible explanation for our observations. In confined one- 
(linionsional ^•ouwtrirs mmk* authors him- pmlicUjd that 



liquid helium might behave as a Luttinger liquid 27, 
For example, in the context of liquid helium-filled car- 
bon nanotubes, Del Maestro et al.[28j have recently used 



FIG. 4. (color online) Values of the fit parameters A and b as 
a function of temperature as determined from data including 
that in Figure [3] (squares), obtained from a stable sample 
studied sequentially over several days; different symbols are 
for data separated by helium transfers to the apparatus. Here 
data that results from | AT | values in the range 10 - 37.5 mK 
have been averaged at each T. For these data sets b is nearly 
independent of temperature with an average value ~ 0.32. 
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FIG. 5. (color online) F determined at (CI + C2)/2 = 
25.6 — 25.8 bar as a function of T for different values of A/t, 
interpolated from the data shown in Figure 3, for the case of 
an applied | AT | = 27 mK. These and data for other imposed 
AT values can be represented by F ~ — ln(T/r). 



quantum Monte Carlo simulations to show that Luttinger 
liquid-like behavior should be present, with a Luttinger 
parameter that depends on the pore diameter. And 
Boninsegni et al.[29j have predicted that Luttinger-like 
behavior will be present in the cores of screw dislocations 
in solid helium. 
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In one dimension, the macroscopic behavior of bosons 
and fermions is the same [30| . A Luttinger liquid is tra- 
ditionally thought of as a one-dimensional fermionic con- 
ductor in which the relaxation of current is due to the 
back scattering of single fermions dressed with bosonic 
phonon-type modes. The bosonic counterpart of this pic- 
ture is the relaxation of a supercurrent due to quantum 
phase slippages 3l|, l32j . In the low temperature limit in 
the presence of a finite chemical potential difference, a 
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